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We present a study of the pseudoscalar propagator in the mixed chiral regime with valence quark 
masses in the e— regime and sea quark masses in the p— regime. We first show the NNLO pre- 
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diction of this observable in the chiral expansion. In sectors of fixed topology, the correlator has 
a pole in \jm\ that can be matched to the topological zero-mode contributions to the correlator. 
I We compute the residue of this pole in a Nf = 2 mixed-action simulation and compare the results 

with the prediction in Chiral Perturbation Theory (ChPT). 

> 

X 

S3 



The 30th International Symposium on Lattice Field Theory 
June 24 - 29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http ://pos.sissa.it/ 



Nf = 2 chiral dynamics in the mixed chiral regime 



Gregory Vulvert 



1. Introduction 

Lattice QCD allows to extract physical results from unphysical conditions such as finite vol- 
ume or unphysical quark masses. In particular the up and down quark masses are typically heavier 
in lattice simulations that in nature. Nevertheless chiral extrapolations allow to extract physical 
results. A more extreme generalization is the partially-quenched setup where quark masses are 
different in the valence and sea sectors. 

It is also possible to use different discretizations in the valence and sea sectors: these are 
the mixed-action formulations [1, 2]. One such setup is to consider overlap valence quarks on 
improved-Wilson sea, so that an exact valence quark chiral symmetry can be preserved at a moder- 
ate cost. This symmetry can be very useful to deal with renormalization of four-fermion operators 
or to access other kinematical regimes such as the e— regime. 

In this work we present the result for the pseudoscalar propagator at NNLO in the mixed- 
regime of ChPT where some quarks are in the e— regime and some in the /?— regime [3, 4], The 
same computation was done with all quarks in the e— regime in [5]. When valence quarks are in the 
£-regime, there are exact poles in \/rrq when the averages are considered in fixed topological sec- 
tors. The residue of these poles can on the one hand be computed in ChPT in terms of low-energy 
couplings, and on the other computed on the lattice from the wave-functions of the topological 
zero-modes. The latter being potentially more efficient than methods requiring the computation of 
the full correlation function [5]. 

We have computed these observables in an exploratory mixed action simulation with valence 
overlap fermions on a sea of Nf non-perturbatively &{a)— improved Wilson fermions (generated 
within the CLS common effort 1 ). We comment on the comparison of the numerical results and the 
expectation from ChPT. 



2. Topological zero-modes wave functions 

Writing the spectral decomposition of the quark propagator 



l = y v,(*)vf 00 + y v*(*)vf 00 (21) 
xy A- mV V {li + m)V' V ' 



zero— modes non zero— modes 



it can be shown that a two-point function, computed at fixed topology v (indicated by (--^v) 
contains a pole in ml due to exact zero-modes. Its residue is given by: 

lim {m y V) 2 Cf (x) = Tr \r a T b ] srf v (x) + Tr [r fl ]Tr [T b ]J v (x) (2.2) 

where 

<(x-y) = (zv)(x)v i (x)vj(y)v j (y)^ , ^ v (x- v) = - ^£vj(x)v,-(x)£v](v)v y (y)^) . (2.3) 

The sums run over the set of |v| zero modes v,- of the Dirac operator, Dv, = 0, which have definite 
chirality, and are assumed to be normalized so that J dx v](x)vi{x) = V. The left side of eq. (2.3) 

1 https://twiki.cem.ch/twiki/bin/view/CLSAVebHome. 
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Figure 1: Illustration of the two types of quark lines contractions. On the left, the connected contraction, 
corresponding to the amplitude g$ . On the right the disconnected one, corresponding to the amplitude ' . 

corresponds to a "connected" contraction of the quark lines whereas the right part corresponds to 
a "disconnected" one. The figure (1) illustrates these features. It is more efficient to compute the 
residues themselves than the full all-to-all propagator, so it makes sense to use these observables in 
the matching to ChPT. 

3. Pseudoscalar correlator in ChPT in the mixed regime 

We consider QCD with N v valence quarks of mass ra v and N s sea quarks of mass ra s . The 
two-point function of the pseudoscalar density that we want to compute is defined by: 

Tr[T a T b ]cf{t) = jdx (f*(x)F*(Q)) ^ , P°{x) = %T a iy 5 ^, (3.1) 

where T a , a = 1, . . . , Ny — 1 is a traceless generator of SU (N v ) and T° = l v . ChPT can predict the 
behaviour of the pseudoscalar correlator in finite volume V = l?T and as a function of the quark 
masses and of a finite number of low-energy couplings. 

In the mixed regime we consider, the valence quarks are in the e— regime and the sea quarks 
in the /?— regime: 

m v LV < 1 and m s EV > 1. (3.2) 
The power counting is thus the following: 

/? 2 ~ZT 2 ~m s ~^(£ 2 ) and m v ~^(e 4 ). (3.3) 

In the £— regime, the only coupling appearing at NNLO is F. However, in the mixed regime, the 
higher order couplings L4, . . . ,Lg contribute because of the heavier sea quarks in the p— regime. 
The sea quarks therefore behave as decoupling particles: we expect that the correlator (3.1) has 
the same behavior as the one in the quenched theory with e— regime quarks, but with effective 
low-energy couplings which are well-defined functions of those of the Nf = N s theory. Indeed this 
expectaction is satisfied. We refer to [3, 4] for more details on the Feynman rules in the mixed- 
regime. 

More concretely, we consider the volume average of the euclidean time derivative of the 
residues of the correlators: 




lim (m Y V) 2 ± [d 3 xC?(x) = U' v (t), 

m v -K) at J 2 

lim (m v V) 2 y / d 3 x C^(x) = N v A' v (t)+N*A' v (t). 
m v ^o dt J 
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dx£/ v (x) andj3^,(?) = J dx£/ v (x) 



The spectral representation of (2.2) and the definitions sa v (t) - 

implie that a large t : sf' v {t) =A' v (t) and ^(t) = A' v {t) . 

We have computed these quantities at NNLO for any N v and N s , but we will focus here on the 
partially-quenched case (N w = 0), which is the relevant one for the numerical study of this work. 



4. ChPT predictions in the partially- quenched theory 

We use the replica method [6] : the perturbative modes give contributions that are finite in the 
limit N y — > 0, while the contribution of the non-perturbative modes involve zero-mode integrals on 
the U (N V \N V ) supergroup. We checked explicitely that the results obtained are consistent with the 
ones obtained directly from the supersymmetric method [7]. 

We obtain the following structure for the correlator: 



C q + a^G(x,0) + a^F(x) [G{x,0)] l + ^ [F {x)f + ff>G(x,0)F (x 



A2), 



[0 



(2) 



+ 



d 4 y 



i 1] G(x - y, 0)G(y, 0) + y\ z, G(x -y,0)F (y) + tf>F (x-y)F (y) 



(4.1) 



+ j3f [G{x,Mi/2)]+e q d^(x), 
with the following definition of the propagators: 



G (x,M 2 ) 



jp-x 



1 



V ^ 4 p 2 +M 27 



°n,0 



G(x) 

e ip-x 
„4 ' 



1 

V 



IN 



(4) 
';/.() 



n 2 ' 



no n 
T'L 



2n 



(4.2) 



(4.3) 



Except for the first term of the third line, the result has the same functional form as the fully- 
quenched result in the e— regime [5]. However the coefficients are different: they are now functions 
of the low-energy couplings of the Nf = N s theory and of the quark masses, as expected. In partic- 
ular, the singlet couplings ttVq and a appearing explicitly in the quenched case do not show up in 
this case since the singlet is integrated out in the partially-quenched theory. However, both results 
match provided the following identifications are used: 

1 - j| [G (0,M 2 /2) - 8L 4 M 2 ] | (4.4) 

1 - i [N S G (0,M 2 /2) - 8L 5 M 2 ] j (4.5) 
X 2 - 1 



1 eft 

Oeff 
2N r 



in, 



O.eff 



2N r 



1 

Ml 
N 



F 2 



1 



1 

F 2 



N 



G (0,M 2 ) - 16M 2 (N,L 6 +N S L 7 +L 8 ) 



(4.6) 



In [5], a large— N c expansion was combined with the chiral expansion so that \/N c ~ ff{e ), 
and the scaling of the singlet couplings was chosen to be i7Vq/N c ~ M 2 ^(e ) and a/N c ~ 0(e 2 ). 
We do not get the same power counting since eqs. (4.5-4.6) imply that Mq e ^/N c ~ G{ 
tteff/Nc ~ ^(1)- To achieve the same power counting we would need to assume l/N s ~ C[e 
Under this assumption, the present results precisely match with the fully-quenched ones of [5], as 
expected. 



and 

2^ 
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A useful quantity to match lattice QCD results is D v so that 



■^A' v (t) = D vZ + C v z 3 + 0(z 5 ), z = ?-\ and t= ^ : 



(4.7) 



and similarly for D v . The study of the quenched case [5] showed that D v can be robustly extracted 
from a fit of this form to the temporal dependance of the time derivative of the correlators. The 
ChPT prediction for D v is: 



2\v\ 



+ 



I i , Oeff , P 



(F eff L) 2 



-ftp 



-3/2 



|^ + ^ + 3 + 2v 2 -2(v 2 > eff 



1 2\v\ 1\ M 2 2 /l A. 4 _, /M 2 



r 4 (73+4C 4 ) 



v|74(M s 2 s /2) 



(4.8) 



m 2 --F 2 

with (v 2 ) e ff = V — e (as in the Witten-Veneziano relation). We define here p = T/L; and the 
numerical values for j3i , Q and for i = 1 , . . . , 4 are given in table (1). 



p 


1 


2 


Coeff\ 






ft 


0.140461 


0.083601 


Ci 


1 


C2 


-1/24 


?3 


7/5760 


C4 


-31/967680 


7i 


-0.057128 


-0.083291 


72 


-0.001951 


-0.002951 


73 


-0.000066 


-0.000101 


74 


pion sea mass dependent 



Table 1: 

Geometry dependent coefficients. 



Lattice 


Total # 
conf. 


# conf. by topological sector 


v= 1 v=2 v=3 v=4 


D4 


117 


43 18 31 25 


D5 


129 


40 43 23 23 


D6 


153 


85 43 25 



Table 2: 

Data sample. 



5. Numerical results 

We have used N s = 2 CLS &(a)— improved gauge configurations on which we built overlap 
valence fermions in a volume of 24 3 x 48 at /3 = 5.3, which corresponds to a value of the lattice 
spacing a = 0.0649(10) fm[8]. Equivalentely, these values correspond to a lattice size L = 1.56 fm 
for p = 2. We use three lattices D4, D5 and D6 corresponding to K = 0.13620, 0.13625 and 
0.136635 respectively. 

The reader can refer to [9] were topological and spectral observables have been studied on 
the same set on configurations. Thus these are already classified by topological charge and we 
could therefore compute the zero-mode saturated correlators at fixed topological charge. Table (2) 
summarizes the data sample. 
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Lattice 


Fitrange 


V = l v=2 v = 3 v=4 


D4 


24-33 
24-35 
24-37 
24-39 


3.0(0.6) 8.9(0.6) 11.5(0.8) 16.6(1.6) 
3.0(0.5) 8.6(0.5) 11.7(0.7) 17.5(1.3) 
2.9(0.5) 8.6(0.5) 12.0(0.7) 18.5(1.2) 
2.9(0.4) 8.9(0.5) 12.3(0.7) 19.5(1.1) 


D5 


24-33 
24-35 
24-37 
24-39 


2.1(0.9) 6.5(0.9) 11.9(1.3) 17.8(1.4) 
2.7(0.7) 7.1(0.8) 13.5(1.0) 19.0(1.0) 
3.1(0.6) 7.7(0.7) 14.6(0.9) 20.0(0.9) 
3.3(0.6) 8.3(0.6) 15.3(0.9) 21.2(0.9) 


D6 


24-33 
24-35 
24-37 
24-39 


4.4(0.6) 9.8(1.0) 16.3(1.3) 
4.5(0.5) 10.0(0.9) 16.7(1.1) 
4.6(0.5) 10.5(0.8) 17.2(1.0) 
4.7(0.4) 11.1(0.7) 17.8(1.0) 



Table 3: 

D v extracted from the linear fit of A'(t) /L 2 for several time intervals. 



Fig. (2) shows the coefficient D v /\v\ as a function of time for two of our lattices and several 
topological charges. This temporal dependence is very well described by the expression given in 
eq. (4.7). We give some results in table (3) for the case of a linear fit of the z dependence of D v . 
We checked the stability of the fits by varying the fitting time interval and also by including higher 
order corrections (with a term in z 3 ). In any case, the results we obtain are in full agreement. 




5 10 15 20 25 30 35 40 45 5 10 15 20 25 30 35 40 45 



Figure 2: Plots of D v /\v\ vs t for two of our lattices (left : D4 - right : D6) for several topological sectors. 

Unfortunately, the convergence of the chiral expansion is bad at this value of L: for L = 2 fm 
and p = 2, for M ss « 250 MeV, the NLO corrections represent 100% of the LO values as illustrated 
by fig. (3). It is therefore not surprising that fits of D v to the NNLO prediction are not good. We 
expect a significantly better convergence for p = 1 and L > 3 fm, with the sea masses we are 
considering. 

6. Conclusions and outlook 

We have studied the pseudoscalar correlator in a mixed regime where the volume is large com- 
pared to the Compton wave-length of sea-quark pions but short compared to that of the valence- 
quark pions. ChPT can be adapted to this regime, and we have shown the result of a NNLO 



6 



Nf = 2 chiral dynamics in the mixed chiral regime 



Gregory Vulvert 







-0.5 









AT* 


= 250McV " 




Af* 


= 300McV " 



2.5 5 7.5 10 

L [fm] 

Figure 3: {D^ LO /\v\- Df /\v\) / (Df /\v\) vs. L. Black and red lines are the results for M ss = 250 and 
300 MeV while dashed and solid curves correspond respectively to p = 1 and 2. 

computation of this observable. In particular we have focused on the coefficient of the 1 jm^ con- 
tributions that appear when considering fixed-topological sectors. At the fundamental level these 
terms can only arise from the topological zero-mode contributions to the correlator. The numerical 
computation of these quantities relies therefore on the evaluation of the zero-mode wave functions, 
the full propagator is not needed. 

We have evaluated for the first time these quantities in a mixed action simulation with overlap 
valence quarks on Nf = 2 improved Wilson sea fermions, with a fine-grained but relatively small 
L = 1.6 fm lattice. We have found that these observables can be extracted robustly and efficiently. 
However a meaningfull matching to ChPT requires significantly larger volumes L > 3 — 4 fm to 
guarranty that the neglected higher order corrections are small. 

The mixed regime considered here can be useful to extract low-energy couplings of ChPT 
in a different kinematical regime with different systematics and different weight of the different 
low-energy couplings (e.g. those that give contributions proportional to the valence quark mass 
would be suppressed). It might also be useful to simplify other notoriously difficult problems, 
such as clarifying the role of the charm quark mass in the Al = 1/2 rule and extracting the weak 
low-energy constants relevant in K — > nn decay. 
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